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Abstract 

Gravitational waves in isotropic cosmologies were recently studied 
using the gauge-invariant approach of Ellis-Bruni We now con- 
struct the linearised metric perturbations of the background Robertson- 
Walker space-time which reproduce the results obtained in that study. 
The analysis carried out here also facilitates an easy comparison with 
Bardeen. 
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1 Introduction 



In a recent paper jjj the gauge-invariant and covariant approach of Elhs- 
Bruni j2j is used to examine shear-free gravitational waves propagating through 
isotropic cosmologies. In this approach the waves are modelled as small per- 
turbations of the Robertson- Walker space-time. The presence of the waves 
is found to perturb the shear and also more notably to introduce anisotropic 
stress into the universe. Other basic gauge-invariant quantities, for example 
the vorticity and energy flow, remain unchanged by the presence of gravita- 
tional radiation. 

Our purpose here is to construct the metric perturbations of the Robertson- 
Walker space-time which give rise to the perturbations of the anisotropic 
stress and shear found in [Ij. The difficulty is that we wish to derive gauge- 
invariant perturbations and there is no way a priori to identify which terms 
in the perturbed metric are pure gauge terms without carrying out a lengthy 
calculation. In the process of studying the perturbed metric we identify the 
gauge terms and without loss of generahty we then put these terms equal to 
zero. 

The paper is organised as follows: In Section 2 we introduce the nota- 
tion used and give some important equations. The unperturbed Robertson- 
Walker space-time is described in Section 3. In Section 4 we summarise the 
results of the gauge-invariant and covariant study of gravitational radiation 
carried out in [1 . The perturbed metric is introduced in Section 5. Also in 
this section and Section 6 we demonstrate how the perturbed metric leads 
to the required gauge-invariant perturbations of the shear and anisotropic 
stress. The Ricci tensor components of the metric are listed in Appendix A 
and in Appendix B we briefly outline the calculation involved in identifying 
those variables which are responsible for the presence of gauge terms. The 
paper ends with a discussion in which our results are compared with those 
of Bardeen |2j . 

2 Notation and Basic Equations 

Throughout this paper we use the notation and sign conventions of 
We are concerned with a four dimensional space-time manifold with metric 
tensor components Qab in a local coordinate system {a;"} and a preferred 
congruence of world-lines tangent to a time-like vector field with components 
u"" and u"'Ua = —1. With respect to this 4-velocity field the symmetric 
energy-momentum-stress tensor T'*'' can be decomposed as 

T"'' = f^u'^u'' + ph'''' + q'^u'' + q^u" + 71"'' , (2.1) 
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where 

h''^ ^ g"^ + u'' u\ (2.2) 

is the projection tensor and 

q^Ua^O, n'^'ua^O, n\^0, (2.3) 

with Here fx is the matter energy density measured by the observer 

with 4-velocity w", p is the isotropic pressure, g'" is the energy flow relative 
to u"' (for example heat flow) and tt"^ is the trace-free anisotropic stress (due 
to processes such as viscosity). 

We indicate covariant difi^erentiation with a semicolon, partial differentia- 
tion by a comma and covariant differentiation in the direction of u"" by a dot. 
Also as usual square brackets denote skew-symmetrization, round brackets 
denote symmetrization and a definition is indicated by a colon followed by 
an equality sign. Thus the 4-acceleration of the time-like congruence is 

= (2.4) 

and Ua;b can be decomposed into 

Ua;b=^^ab + (^ab+^Ohab-UaUb , (2.5) 

where 

i^ab-^Ula-b]+U[aUb], (2.6) 

is the vorticity tensor of the congruence, 

Crab-^U(^a;b)+U(aUb)-^Ohab, (2.7) 

is the shear tensor of the congruence and 

^:=<a, (2.8) 

is the expansion (or contraction) of the congruence. 
We shall make use of the Ricci identities 

Ua;dc — Ua-cd = RabcdU^ 1 (2-9) 

where Rahcd is the Riemann curvature tensor but for the problem at hand 
the key equations are Einstein's field equations 

Rab-\9abR = Tab . (2.10) 
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Here Rat '■= Ra^bc ^^re the components of the Ricci tensor, R := R'^c is the 
Ricci scalar and we have absorbed the couphng constant into the energy- 
momentum-stress tensor. Noting that R = —T{:= T°'a) and using Eq. 
()2.H) the field equations can be decomposed into: 

RabWv!' = i(/i + 3p) , 

RabW'h'^ = -g, , (2.11) 
Rab Kh!'^ = 2 ~ ^'"^ + '^'"^ ■ 
It is in this form that we shall use Eq. ()2.1()|1 in later sections. 

3 The Background Space— Time 

We choose as the unperturbed (background) space-time a Robertson- 
Walker space-time with line-element 

,s^ = n^,,M^ll±ldEl±(^.,,^, (3^) 

(l + |r2) 

where R{t) is the scale factor, = {x-^^ + (x^)^ + (x^)^ and k = 0, ±1 
is the Gaussian curvature of the space-like hypersurfaces t = const. The 
world-lines of the fiuid particles are the integral curves of the vector field 
u°'d/dx"' = d/dt (thus = 6^ since we shall label the coordinates = 
y ,x'^ = z ,x^ = X ,x'^ = t)- The background energy-momentum-stress 
tensor is Eq. ()2.H) specialized to a perfect fiuid (by putting q°- = = tt"^) 
with proper-density 

and isotropic pressure 

P = 

We find it convenient to put the line-element given above in the following 
forms: 

rfs^ = R\t) {dx'^ + f {dy^ + dz'^)} - dt^ , (3.4) 

with po = 1 + {.K/A){y'^ + z^), K = const, / = f{x). We identify three 
distinct cases: 
Case 1: 



R^ k_ 



'R^ R W 



(3.2) 



(3.3) 
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If k = +1 then K = +1 and f{x) = sinx. 

Noting that the transformation x 7c/2 — x does not affect the form of the 
hne-element ()3.4|1 we see that in this case f{x) could equivalently be written 
f{x) = cosx 

Case 2: 

K = and f{x) = 1 , 
If = then i or 

K = +1 and f{x) = x . 



Case 3: 

K = —1 and f{x) = coshx 
or 



li k = —1 then < 



K = and /(x) = | 



2 

or 

i^T = +1 and f{x) = sinhx . 
The form of the hne-element ()3.4j) is also invariant under the transformation 
X —X so when if = in case 3 we could instead write /(x) = | e~^. For a 
detailed explanation why these cases arise see for example Eqs. (5.3)-(5.19) 
in p. In space-times with line-elements ()3.4jl the hypersurfaces 

0(a;") ■=x- T{t) = const , (3.5) 

with dT/dt = are null hypersurfaces. The expansion of the null geodesic 
generators of these surfaces is 

1 f'R 

where /' = df /dx, R = dR/dt. Using ()3.5|) we can show that 

20,a;6 = (t),b + 6 0,a + ^.d'*^ Qab , (3.7) 

where 

f' 

^a = -j(t>,a + R<P/Ua. (3.8) 

It follows from Eq. ()3.7|1 that (j)^a is shear-free [Hj. 

Finally in this section we note that for convenience we have used the 
same coordinate labels {?/, t} for all the special cases included in ()3.4p . 
Clearly the ranges of some of these coordinates will vary from case to case 
and within cases 2 and 3. For example, in case 2 x G (— oo , +oo) if i^' = 
but X e [0, +cxo) and is a radial polar coordinate ii K = +1. The shear-free 
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null hypersurfaces ()3.5|) will also be different in the different cases. This can 
be seen by examining the intersections of these null hypersurfaces with the 
space-like hypersurfaces t = const. 
Case 1 

The intersection is a 2-sphere. 
Case 2 

If i^' = +1 the intersection is a 2-sphere and if i^' = the intersection is a 
2-plane. Thus it is obvious that Eq. ()H.5j) describes two different families of 
shear-free null hypersurfaces that can occur in an open, spatially flat uni- 
verse. 

Case 3 

In this case the intersection of with the t = const hypersurfaces is al- 
ways a 2-space of constant curvature. The curvature of this 2-space is given 
by K which takes values 0, ±1. So we have three different families of shear- 
free null hypersurfaces in a A; = —1 universe. We refer the reader to jH] for a 
geometrical explanation for the existence of these subcases. 

4 Gauge— Invariant and Covariant Approach 
to Gravitational Waves 

In a recent paper P we used the gauge-invariant and covariant ap- 
proach of Ellis-Bruni |^ to construct gravitational wave perturbations of 
the Robertson-Walker space-times described in the previous section. This 
involves working in a general local coordinate system with gauge-invariant 
small quantities which by their nature vanish in the background, rather than 
small perturbations of the background metric. For isotropic space-times the 
Ellis-Bruni variables are aab, u", ujab, Xa = hlfi^b, Ya = Kp^b, Za = h^O^b, 
T^ab, Qa and the "electric" and "magnetic" parts of the Weyl tensor, with 
components Cabcd, given respectively by 

Eab = Capb, uP , Hab = '' Capb, u" . (4.1) 

Here *Capbq = \ Vap^^ Crsbq IS the dual of the Weyl tensor (the left and right 
duals being equal), r]abcd = y/^ ^abcd where g = det(gab) and eabcd is the Levi- 
Civita permutation symbol. However we found that it is tensor quantities 
that describe gravitational wave perturbations. Thus for this problem the 
important Ellis-Bruni variables are aab, T^ab, Eab, Hab and we can set all other 
gauge-invariant variables equal to zero. The equations satisfied by these 
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variables are obtained by projections in the direction u"" and orthogonal to 
u"" of the Ricci identities, the equations of motion and the energy conservation 
equation contained in T"'^-b = and the Bianchi identities written in the form 

C^bcd^^ = R4a;b] _ 1 gc[a j^,b] _ f^^ ^) 

To keep this section to a reasonable length we shall not list all of the equations 
(they are given in Eqs. (2.14)-(2.25) in P). We note here that from the 
projections of the Ricci identities (after putting -u" = = Uab) we find 

12 2 

Eab = -TTab + -<y'^ hab - - CFab " C^af b " H hf, dab , (4.3) 



and 



Hab = -Khla^t''''r]s)fgcuf . (4.4) 



Thus these variables are derived from nab and aab- 

We now assume that the perturbed shear and anisotropic stress have the 
following form: 

CTab = SabF{(f)) , TT^i, = UabF{(f)) , (4.5) 

where F is an arbitrary real-valued function of its argument (/)(a;°). We 
emphasise that at this point 0(a;") is arbitrary and not that defined in Eq. 
()3.5p . This idea of introducing arbitrary functions into solutions of Einstein's 
equations describing gravitational waves goes back to work by Trautman 
[Zj and the above form for the gauge-invariant variables was introduced by 
Hogan and Ellis jH]. Substituting ()4.5|) into the linearised versions of the 
equations satisfied by these variables and noting that Sab and Uab are trace- 
free and orthogonal to Ua with respect to the background metric we find that 

in 

9''' (t^,a 0,6 = , s'^' = , n-^' = , (4.6) 
with Qab here the background metric, and 

s'^V = , = , (4.7) 

where for clarity we have used a stroke to denote covariant differentiation with 
respect to the background metric. We also discover (see [T]) the following 
wave equation for Sab 

^a6|d ^ _'^Q^ab _f}_Q^^Q2\ ^ab + (p _ 1 ^) g'^' = -fl"^ - ^ ^ H"^' , (4.8) 
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and a propagation equation for Sab along the null geodesies tangent to 0'"^, 
namely, 

4 + [l - ^ ^ 0) = -i n,, , (4.9) 

where s[f, := Stb\d4>''^ and = (f),au"'- The internal consistencies of these 
equations were checked in pp. The "electric" and "magnetic" parts of the 
Weyl tensor are now given by [T] 

Eab= {^^llab-Sab-'^eSab^ F - <P S^b F' , (4.10) 

and 

Hab = -S(a^'' Vb)fpc uf F - S(/ rib)fpc <P'^ F' , (4.11) 

where F' = dF/ dcj). These equations are easily checked by substituting ()4.5|1 
into (jOI) and (jOjl . 

We wish to construct pure gravitational wave perturbations i.e. having 
pure type N perturbed Weyl tensor in the Petrov classification. It is shown in 
jl] that on account of ()4.6|1 the F -parts of Eab and Hab above are type N with 
degenerate principal null direction 0'"^. Then if we also require the F parts of 
Eab and Hab to be type N the perturbations we have constructed describe pure 
gravitational waves with propagation direction cf)'^ in the Robertson-Walker 
background and the histories of the wave-fronts are the null hypersurfaces 
0(x'*) = const. Making use of the following null tetrad, ka = —4>~^ (p,a, 
la = Ua — \ ka and m^, fha a complex covariant vector field and its complex 
conjugate chosen so they are null (m" iria = = fff- fha)., are orthogonal to fc" 
and 1°" and satisfy m° fha = 1 we find that a simple way to ensure the F-parts 
of Eab and Hab are type N is to require the null hypersurfaces 0(x") = const 
to satisfy (see P) 

0,fe|em'r = O, (4.12) 

and 

0,„|fem"m^ = . (4.13) 

To exhibit explicit examples we specialise to the case = x — T(t) with 
T{t) introduced in ()3.5|1 . Then the null tetrad described above is given by 
the 1-forms 

ka dx"" = Rdx — dt , la dx"' = — {Rdx + dt) , 
rriadx^ = -^RpQ^f{dy + idz), (4.14) 

and it is straightforward to check that Eqs. ()4.12j) and ()4.13j) are satisfied. 
Since s"* and 11"'' are trace-free and orthogonal to and 0'", they each have 
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only two independent components. These components are = —s^^ = 
a{y,z,x,t), = ^21 ^ (3{y^z,x,t) and = -H" = A{y,z,x,t), li^^ = 
iP^ = B{y, z,x,t) where we have labelled the coordinates x^ = ?/ = 
z ,x^ = X ,x'^ = t. Now we can write 

s^^ = sm^'m^ + STh''m\ (4.15) 

with 

s = -i?Vf + , (4.16) 

and 

n'*'' = nm"m^ + nm"m^ (4.17) 

with 

n = -R^pff{A + iB) . (4.18) 

It follows from ()4.7j) that a, (3 and A, B must satisfy the Cauchy-Riemann 
equations 

^(Po'«) - = , (4.19) 

^(Po'/3) + ^(Po'«) = 0. (4.20) 

and 

Yy^Po'A) - §;iPo'B) = , (4.21) 

§^(Po'B) + §-^ip,'A) = 0. (4.22) 

If we define Q = Pq^ p B?{a + i (3) and note that f = f{x), R = R{t) and a, 
(3 satisfy Eqs. ()4.19j) and ()4.2()j) then Q is an analytic function of C, := y + iz. 
We can now rewrite Eq. ()4.1(ij) as 

s = -R~Vj-'g{Cx,t) . (4.23) 

From the propagation equation ()4.9j) we find 

n = -2 R'yj-\Dg + Rg) , (4.24) 

where D is given hj D = d/dx + Rd/dt = d/dx+d/dT and the dot indicates 
differentiation with respect to t. As a consequence of this and ()4.18p A + iB 
is analytic in ( and so Eqs. 1)4.211) and ()4.22p are automatically satisfied. 
Replacing s""^ by Eqs. and and by Eqs. pTfjl and ^TI^ 

the wave equation ()4.8|1 simplifies to 

D^g + kg = 0, (4.25) 
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with /c = 0, ±1 labelling the Robertson-Walker backgrounds with line-elements 
of the form ()3.4|) . The solutions of these three differential equations are: 
for k = 0, 

g{C, X, t) = a{C, x-T){x + T) + b{C, x-T) , (4.26) 

for k = +1, 

GiC, X, t) = a(C, x-T) sin + b{C, x-T) cos , (4.27) 

and for k = —1, 

e;(C, X, t) = a(C, x-T)sinh(^±^)+6(C, x-T) cosh , (4.28) 

where in each case a{(, x — T), b{(, x — T) are arbitrary functions. Using the 
identity x + T = 2x — (x — T), (and some simple trigonometric and hyperbolic 
relations) we can rewrite ()4.26j) in the form 

GiC, X, t) = h (C, x-T)+ xhiC, X-T), (4.29) 

with hi, h2 arbitrary, ()4.27j) as 

X, t) = /i3(C, X — T) sinx + /i4(C, x — T) cosx , (4.30) 

with h^, h^^ arbitrary and ()4.28|) as 

^(C, X, t) = h^i^C, x — T) sinhx + hQ{(, x — T) coshx , (4-31) 

with h^, arbitrary. In addition 1)4.311) can be put in the form 

GiC, X, t) = hriC, X - T)e^ + h^{C, x - T)e-^ . (4.32) 

When these results are derived from metric perturbations in Section 5 below 
the expressions ()4.29p - ()4.32p will be more useful for comparison purposes 
than the equivalent expressions ()4.26jl - ()4.28|) . 

The "electric" and "magnetic" parts of the Weyl tensor (Eqs. (14.10)1 and 
1)4.11)1 respectively) are now calculated and we find that they can be written 
compactly as P 

^ab j^iH"^ = -2 R-^plf-^-^ig F) m" . (4.33) 

ox 

Here Q is given by Eqs. (jOHIl-dOHI) (or equivalently ^^-^^) and 
F = F{x- T) so that F' = dF/dx, po = I + {K/4:){y'^ + z"^), f = f{x) 
described in the previous section and R{t) is the scale factor. It follows from 
Eqs. 01311, (jug), dmZjl, ()0!?|l and ^^Mi that to find aat, T^ah from Sah 
and Ilab we simply replace Q hj Q F. This does not affect Eqs. ()4.24)1 and 
1)4.25)1 since DF = 0. Conversely with F = F{x — T) and Q given by Eqs. 
dOnil-dOHl) F can be absorbed into Q. 
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5 The Perturbed Metric 

We now exhibit a line-element which (i) can be viewed as a perturbation 
of the space-time line-element ()H.4|1 and (ii) produces the same explicit per- 
turbations described in the gauge-invariant formalism of the previous section. 
We first introduce a pair of null coordinates, 

« = -^(x-T(t)) , v = l=(x + T{t)) , (5.1) 

with T{t) introduced after Writing 

R{t{T)) = n{T) = n{v - u) , (5.2) 

the line-element ()3.4|1 written in terms of u and v reads 

ds'^ = fidy"^ + dz'^) + 2n^dudv , (5.3) 

where now / = f{u + v). The coordinates y, z, u, v are such that the surfaces 
u = const, V = const are two families of intersecting null hypersurfaces. The 
general form of line-element in a coordinate system based upon two families 
of intersecting null hypersurfaces is given in For our purposes we write 
this as 

ds"^ = 6^ hj^B (dx^ + af du + dv) {dx^ + du + dv) + 2cdudv , (5.4) 

where A, B take values (1,2), {hAB{y, z, u, f )) is a unimodular 2x2 symmetric 
matrix, (x\x^) = {y,z) and af, a^, b, c are six functions of y,z,u,v. It is 
convenient to use the following parametrisation JU] of {hAB)'- 

, /e^" cosh 2/? sinh2/3 \ . . 

^^^^)-[ smh2p cosh2/? j ■ ^^'^^ 

Here a, j3 are taken to be small of first order. With (/iab) given by Eq. ()5.5|) 
it is easy to check that, working to first order, Eq. ()5.4|) can be written 

ds"^ = b'^[{l + a)dy + (3dz + {al{l + a) + all3}du + {al{l + a) + al(3}dvf 
+b'^[(3dy + (1 - a)dz + {a} /3 + a^(l - a)} du + {all3 + al{l - a)}dvf 
+2cdu dv . (5.6) 

The background space-time is obtained from this by putting 

af = , a^ = 0, b = pQ^nf , c = , a = , p = . 

(5.7) 
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For the perturbed space-time that we require we find that 6, c retain their 
background values, and we can put = = a^. These latter quantities 
actually play the role of gauge terms (see Section 7 below and Appendix B 
for an illustration of this). Every shear-free system of gravitational waves 
involves an arbitrary analytic function pTj and we now have two real func- 
tions a, j3 available to provide the real and imaginary parts of this analytic 
function. Also we find that the 4-velocity u"", the isotropic pressure p and 
the matter-energy density yU take their background values (these are given 
in Section 3). 

To demonstrate that this space-time does indeed describe the perturba- 
tions of Section 4 we shall work on the tetrad given via the 1-forms 

= p^^ fVl{{l + a)dy + pdz} , 

9^ = p,^ fn{l3dy+{l-a)dz} , 

9^ = ndu , 

9^ = VLdv , (5.8) 

with po = 1 + {K/4:){y'^ + z^) as in ()3.4|1 . We note that with respect to this 
tetrad the line-element is now 

ds^ = {9^f + {9^ + 29^9^ = gab , (5.9) 

thus defining the tetrad components Qab of the metric tensor. The tetrad 
components of the matter 4-velocity are given via the 1-form 

Ua^'^ = ^(^'-^') . (5.10) 

Since we wish to reproduce the linear perturbations of the previous section 
we shall discard any terms which are second order or smaller in a and (3. Our 
first step is to calculate the Ricci rotation coefficients and the Ricci tensor 
components. This results in a lengthy list of equations which for convenience 
we give in Appendix A. We now use the Ricci tensor components and fj5.10|) 
in the field equations given by Eqs. ()2.1H) . Noting that 

Q' = ^RR, Q" = -R'^ R+-RR'^ , (5.11) 

where the prime denotes differentiation with respect to v and using Eqs. 
()3.2j) and ()3.3j] it is easily checked that the first of Eqs. ()2.1H1 is identically 
satisfied. The second equation in ()2.11|) yields 

3 

qi = ^ f-'n-'{ip^'a)yu + {Po'P)zu - {pfa)y, - {pfP).v} ,(5.12) 
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52 = ^ f-'n-^{{Po'a),, - {pfl3)y, - {pfa),u + {pfP)yu} ,(5.13) 

gs = 0, (5.14) 
g4 = . (5.15) 

The subscripts u, v here indicate partial differentiation with respect to 
these variables. We recall that in the covariant approach we found that 
Qa = 0. With a, P chosen so they satisfy the Cauchy-Riemann equations in 
the form: 

ip^^a)y+ip^^f3), = , (5.16) 

{Po'a),-{p^^f3)y = , (5.17) 

it follows that = as required. For later use we note that as a result of 
these equations Po ^(a — i P) is analytic in ( = y + iz. With Qa = the last 
of Eqs. ()2.11|) can be rewritten as 

Rab = 1^ UaUb +phab + TTab-^i^P- f^)9ab ■ (5.18) 

With the Ricci tensor components given by (A.11)-(A.20) it follows from 
this and Eqs. (Q, (Q, (EH)- (IHTT|l . (irTH|l . (IKTfl) that TCab = except 
for TTii, 7122 and 7Ci2 with 

7iu = V2R-^R{a,-au)-2R-^ f-^ / {a, + au)- 2 R~^a,u , (5.19) 
7T22 = -V2R'^R{a,-au)+2R-^ f-' f'{a, + au)+2R-^a,ui^.20) 
71,2 = V2R-^R{P,-Pu)-2R-'f-'f'{P, + Pu)-2R-'p,^. (5.21) 

We have made use of 

2f" = -kf, 2{f'f + kp = K, (5.22) 

to simplify these equations. We note that the prime here denotes differenti- 
ation with respect to f = (1/a/2)(x + T). Similar equations to these appear 
in (eq. (5.41)). In ^ the prime indicates differentiation with respect to 
X and hence the factors of 2 in ()5.22|1 do not appear there. 

Now in terms of the background null tetrad described by Eq. ()4.14p we 
can write the coordinate components of the (small) anisotropic stress tensor 
as 

TT*-?' = 7f m-' + TT rh^ . (5.23) 

Using Eqs. ()4.14j) and ()5.8j) (with a = (3 = Q since vTah is a first order 
quantity) we find 

* = o(^ii ~ ^22) - «7ri2 . (5.24) 
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Substituting from Eqs. ()5.19|) -f5.21) above yields 

7f = -2pl [pq^ f{a - i(3)uv + Po"^ f' {av + Oiu - i{Pv + Pu)}] 

-V2R'^ R{au -a,- i{(5u - (^v)} ■ (5.25) 

In the previous section we worked with J\ah and 11 which we expressed in 
terms of an analytic function Q. But i^oh = ^ab F and as indicated following 
Eq. ()4.33p F can be absorbed into Q. Hence, in order to make contact with 
the gauge-invariant description yf here must satisfy the same equation as 11 
and thus we require 

Ti = -2plf-^R-^{Dg + Rg} , (5.26) 
with D = d/dx + Rd/dt = \^d/dv for some analytic function Q. Taking 

Q = -^Po"^ f{<yu- ay -i{f3u- M} , (5.27) 

we find that it is indeed possible to write vf in this form provided we choose 
a, P to satisfy the following: 

If /' = then = ,(5yy = ] (5.28) 

if /V then = , = . (5.29) 

We note that the first of these conditions corresponds to the case k = 0, 
K = described following Eq. ()3.4j) . We now assume that these conditions 
hold. As a consequence of these and Eq. ()5.22j) it immediately follows that Q 
given by ()5.27|) satisfies the wave equation (j4.25|) . Also noting that f = f{x) 
and using the Cauchy-Riemann equations ()5.16|) - ()5.17|) we see that as before 
Q is an analytic function of ( = y + iz. 

We now turn our attention to the shear. In a similar fashion to the 
anisotropic stress the coordinate components of the (small) shear tensor can 
be written in the form 

CTjj = Sp rrii rrij + Sp rhi fhj , (5.30) 

where m,, rhj are given by Eq. ()4.14|) and in terms of the Ricci rotation 
coefficients 

Sp = ^^{(Ti4i-T242 -Ti3i + T232)+«(Ti32 + T231-Ti42-T24l)} . (5.31) 

Evaluating this using the Ricci rotation coefficients given in Appendix A we 
find 

= ^ R^'{{a - iP)y -{a- iP)u} = -Po ' R-' Q , (5.32) 
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with Q as before. Taking into account that we can absorb F into Q and that 
Sp = s F (with s defined by Eq. ()4.15|) ) we see that the perturbations we 
have produced here also satisfy Eq. ()4.23|) . Thus we have shown that the 
perturbations described by the metric ()5.6j) take the same form as those found 
by the covariant approach. In the next section we shall illustrate that they 
also satisfy the wave equation ()4.8|) and the propagation equation ()4.9|) . 

For the remainder of this section we compare the explicit Q found here 
with the solutions of the wave equation found in ^ and listed in Eqs. ()4.26|1 - 
()4.28|1 (or equivalently ()4.29|l - ()4.H2|l ). We first examine the case when k = 0. 
There are two subcases to consider here (i) K = and f{x) = 1, (ii) K = +1 
and /(x) = X. When K = 0, pq = 1 and Eq. ()5.27|) reads 

Q = ^{{a - iP)u - {a - tP),} . (5.33) 

Since /' = in this case we have a^^ = = (3yy. Thus in addition to a — ip 
being analytic in ( this complex-valued function is also linear in v. Hence 
we can write 

g{C, X, t) = ai(C, x-T){x + T) + a^iC, x - T) , (5.34) 

where ai, 02 are arbitrary (analytic) functions of their arguments. When 
K = +1, f{x) = and from ()5.29|1 we have ay = = [3^. Therefore the 
function Po ^(a — i/5) is analytic in ( and independent of f = (x + T)/ v^, i.e. 
it depends only on ( and u = {x — T)/y/2, and we can write ()5.27|) in the 
form 

Q = xa3{C, X — T) . (5.35) 

with 03 an arbitrary analytic function. Using the identity x+T = 2x — {x—T) 
as in Section 4 we can rewrite ()5.34j) in the form ()4.29|1 . Then ()5.35j) is the 
special case of ()4.29|1 corresponding to hi{(,x — T) = 0. Thus in the case 
k = there are two independent expressions for Q{(,x,t) which are given in 
the form of a superposition in ()4.29|) . This arises because ()4.29|) is obtained 
by solving the linear wave equation ()4.25|1 with A; = and in general this 
equation is insensitive to the allowable values of i^' = 0, ±1. 

We now look at the solution when k = +1. There is only one case to 
consider here, K = +1 with /(x) = sinx or equivalently (see Section 3) 
f{x) = cosx. Again we have Po^(a — ip) analytic in ( and independent of v 
so ()5.27p can now be written 

Q = ci4(C? X — T) sinx or Q = a^i^^, x — T) cosx , (5.36) 

where 04, 05 are arbitrary functions. The two equations in ()5.3(ij) are equiv- 
alent to dOOl). 



15 



Finally when k = —1 there are three subcases to look at corresponding 
to i^' = 0, ±1. In all cases / (x) 7^ and so we have Po ^(a — independent 
of V and Q has the form 

g = f{x)ae{C,x-T) , (5.37) 

where is an arbitrary analytic function. When K = we have f{x) = 
or equivalently /(x) = |e~^ and so in this case ()5.37|1 agrees with ()4.32j) . 
When K = +1, f{x) = sinhx and now ()5.37|1 agrees with 1)4.311) when 
Hq^Cx — T) = 0. When K = —1, /(x) = coshx and ()5.37p agrees with 
()4.31|) when h^{(, x — T) = 0. This case k = —1 is a good illustration of the 
insensitivity of the expressions ()4.3H) and ()4.32)1 to the values of K. 

Thus all of the solutions found here are identical to the solutions found 
using the gauge-invariant and covariant approach to perturbations in PJ. 



6 Properties of the Shear and Anisotropic 
Stress 

In the previous section we exhibited a perturbation of the Robert son- Walker 
background line-element ()3.4|) that produced perturbations in the shear and 
anisotropic stress tensors which satisfied some of the equations found us- 
ing the gauge-invariant and covariant approach of Section 4. We now show 
that these perturbations satisfy the remaining equations, namely that the 
anisotropic stress and shear tensors are trace-free, orthogonal to u"", divergence- 
free with respect to the background metric and also satisfy the wave equation 
()4.8p and propagation equation ()4.9|) . To do this we shall, in this section, 
work in coordinate components [in the local coordinates y, z, x, t] instead of 
the tetrad components we have used up to this point. In terms of this local 
coordinate system we can write the line-element ()5.6|) (with = = 0, 
b = Po^fl f,c= in the form 



ds'^ = gab dx"" dx^ + 2'yab dx"' dx^ := gab dx"" dx^ 



(6.1: 



where cjab = diagjpo ^ ^^,]?o ^ ~1} ^^e metric of the back- 
ground space-time and 



(a (3 

l3 —a 



V 



\ 





0/ 



(6.2) 
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is the perturbation. Clearly '-jab is trace-free and orthogonal to u"" = 6^ (with 
QabU'^u^ = — !)• The non-vanishing Christoffel symbols of the background 
metric tensor given via the line-element ()3.4|1 are: 

'-12 ~ 11 ~ 22 ~ 2 ^0 ' 

■pi _ -^2 _ -pS _ Q-1 Q 
■■-14 — ■'-24 — -"-34 — ^ ' ^ ! 

fn = Tt2=Po'f^^t, (6.3) 

Til = ^22 ~ ~Po f fx 1 
-^13 = 1^23 = / /x , 

= ^ . 

We have used the hat here to emphasise that these are background Christoffel 
symbols and we shall continue to use this notation to denote background 
quantities for the remainder of this section. Using these and the Cauchy- 
Riemann equations ()5.1fij) - ()5.17j) it is a simple exercise to show that •jab 
defined above is divergence-free. 

In order to show that it"'' is also divergence-free we first write it in terms of 
lab- We define the perturbation of the Christoffel symbols to be 6 T^^ := T^^ — 
r^^. Noting that = g"'^ — (here 7*^* = g"''^ g''-^ 'jcf and we are neglecting 
second order small quantities ) it is easily derived from the definition of the 
Christoffel symbols that 

SrZ=\{l^\, + r,\b-9''^lM\f) , (6.4) 

where as usual the stroke indicates differentiation with respect to the back- 
ground metric. Now 7"* is divergence-free and thus we can see from this 
equation that ^T^^ = 0. In general the components of the Ricci tensor of a 
perturbed metric can be written in the form 

Rbd = Rbd + iSTl)^a'{5Tl)\,. (6.5) 

For the problem at hand we have 

Rbd = fJ'UbUd + pfiM-'^iSp- fi)gbd ■ (6.6) 

Substituting for Rm and Rbd from Eqs. ()5.18j) and ()().() |1 respectively in Eq. 
dnUl) yields 

{SKd)\a = '^bd + -{f^-phbd ■ (6.7) 
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Taking the divergence of Eq. ()6.4|) and using this equation we arrive at 

1b\da + ld\ba - 9^hbd\fa = SvT^ + (/i - p) Jbd ■ (6.8) 

Next making use of the Ricci identities 

labldc - 1ab\cd = Rafcd lb " ^fbcd l[ , (6.9) 

and recalhng that 7"'' is divergenceless and orthogonal to u"" we find (since 
Cabcd = in the background) 

3 11 1 

7a\dc = 1^11 Red +- Id Rac - -QadRfcl^' - -Rlad . (6.10) 

With R = /i — 3p and Rab given by Eq. ()5.18|) this equation allows us to 
write (since we are concerned here with first-order terms only) 

lb\da + ld\ba = - P) Ibd , (6.11) 

and hence Eq. ()6.8|) now becomes 

2 

9''hbd\fa - ^^'lbd = -2 7rfcd . (6.12) 

It is easy to see from this that T^ab is trace-free and orthogonal to m". Starting 
with this equation we shall now prove that i^ab is indeed divergence-free. 
That this is necessary to fully make contact with the gauge-invariant and 
covariant approach of Section 4 follows from the fact that in this case we 
wrote Tiab = T^ab F with = 0, F = F(x - T) and = except for 
IP^ and and therefore Il°'''\b = is equivalent to vr'^^it = in this case. 
First making use of the Ricci identities for a tensor of type (3,0), Eq. ()2.1|) 
and Eq. ()2.10|) we can write 

r'^'\db = (7"^%)"^ +{1^'-Ip) ^""'^ = ■ (6.13) 

Also since for the perturbed space-time we are considering here the matter 
density fi retains its background value we have h'^^fi^ = from which it 
follows that 

fi^b = -fiUb . (6.14) 

As a consequence of these last two equations we find, on taking the divergence 
of Eq. dnill, that 

n'''ib = 0, (6.15) 
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as required. 

We shall now examine the properties of the shear aab- As with the 
anisotropic stress above it is necessary to express this in terms of 7'*''. This 
is easily done using the definition of the covariant derivative of Ua- 

Ua;b ■■= -Kb Uc = -Kb Kb^c ■ (6.16) 

We remind the reader that the semicolon here indicates covariant differenti- 
ation with respect to the perturbed metric (background plus a small pertur- 
bation) while a stroke denotes covariant differentiation with respect to the 
background metric. In the background Robert son- Walker space-time the 
shear, vorticity and the 4-acceleration all vanish and so Eq. fl2.5|l specialises 
to ^ 

KbUc-=Ua\b = -Ohab, (6.17) 

in this case. Making use of this equation and Eq. ()6.4|) in Eq. ()6.16|) it 
follows, on account of hab = hab + 'Jab (since for the problem at hand is 
unperturbed), that 

Ua-b= -Ohab + -jab ■ (6.18) 

Here and for the remainder of this section a dot indicates covariant differenti- 
ation with respect to the background metric in the direction of u"". Recalling 
that the 4-acceleration is zero in the background Robertson- Walker space- 
time (i.e. Ua\bU^ = 0) it is trivial to see from the latter equation that the 
4-acceleration in the perturbed space-time also vanishes. We also note that 
this equation is symmetric in (a, b) and thus it is clear from Eq. ()2.6|) that, 
as in the covariant approach, the vorticity tensor vanishes in the perturbed 
space-time. Now equating Eqs. ()6.18|1 and ()2.5|l with the 4-acceleration and 
vorticity tensor both zero we arrive at a simple relationship between aab and 
7a6 namely, 

1 

CTab = -iab ■ (6.19) 

Using this and the properties of 'jab it is straightforward to check that aab 
is trace-free and orthogonal to m". However further calculation is necessary 
to show that it is also divergence-free (this is required for similar reasons to 
those given above while discussing the anisotropic stress). First using the 
Ricci identities given in Eq. ()fj.9|l and noting that Cabcd = we calculate 

rV = (7"V)- + \ l"^ Rfc u'-\ Rbf l'^ . (6.20) 

Replacing Rab here by the right-hand side of Eq. ()6.6|) and keeping in mind 
that 7"^% = 0, -i^^Ub = leads to 

rV = ' (6-21) 
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and therefore result of Eq. = 0. 

At this point all that remains to fully make contact with the gauge- 
invariant and covariant description of gravitational wave perturbations out- 
lined in Section 4 is to reconstruct the wave equation ()4.8j) and the propa- 
gation equation ()4.9|1 . This is done as follows: Using the Ricci identities the 
covariant derivative in the direction of of Eq. ()6.12|) can be written as 

-277'^^ = t'^'I'm^^-^^ (p+^/i) 7"' + ^(/i + 3p)r' 

-If^r'-l^^r' ■ (6.22) 



Also with h^fi^b = 0, h^p^b = 0, and h^9^b = we find, again using the Ricci 

2 



identities and Eq. ()(i.l2|l . that 



~^e[^^,Y'-2n'^')-ler'' . (6.23) 
Entering this into Eq. ()6.22|1 and replacing 7"'' by 2 cr"^ we arrive at 



3 '3 V3 9 / V 3 

We have made use of the background values of 6 and fi to write the equation 
in this form. The background value of 6 is 

= -le'-l{f^ + 3p) (6.25) 

which is obtained by specialising Raychaudhuri's equation to the background 
(i.e. putting g^, Hab, Ua, o'ab and Uab all equal to zero) and the background 
value of /i is given by 

= -e ifi + p) . (6.26) 

This is found by specialising to the background the projections along and 
orthogonal to m° of the conservation equation T"*|f, = (see for example 
Eqs. (2.20) and (2.21) in [1 ). Both the wave equation and the propagation 
equation are actually contained in Eq. ()6.24|) . To confirm this we again put 

C^ab = SabF{(t)) , Tlab = IlabF{(f)) , (6.27) 

where F{(f)) is an arbitrary analytic function of its argument (j) = x — T{t). 
In the covariant approach we found s""^ (j)^b = and U"'^ (p^b = 0. This is also 
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true here since (p^b = (0, 0, 1, —R~^) and we have vr'^'' = 0, vr^^ = 0, a^^ = 
and a^^ = 0. In addition since the hypersurfaces 0(a;") = const are null we 
have 4>''^(p,d = 0. Thus we can write 

tt"^ = ri"^F + 0n"^F' , (6.28) 

^ s^^ F + <p s'''' F' , (6.29) 

and 

^abid ^ = s'^^l'^i, F + {2 s'^"!'^ + s'^' 0'^|rf) F' , (6.30) 

where F' = dF/d(f). Substituting these expressions for a"* and vr"^ into Eq. 
()6.24|) and equating the F and F parts separately yields the required wave 
equation ()4.8|) and propagation equation ()4.9|) . 

7 Discussion 

We have shown in Sections 5 and 6 that the perturbations of the background 
Robertson-Walker space-time derived here from metric perturbations are 
exactly the same as those obtained using the covariant approach. Thus the 
metric ()5.6j) with = 0, = and a, (3 chosen to satisfy the Cauchy- 
Riemann equations ()5.16|) - ()5.17|) is indeed that which we set out to find. We 
mentioned earlier that the functions a^, play the role of gauge terms. That 
this is true is seen by repeating the calculation of Q with 7^ 0, 7^ 0. To 
save repetition here this calculation is outlined briefly in Appendix B. The 
result is that , do not appear in the required analytic function Q i.e 
that which satisfies Eq. ()5.26|) . Thus since all gauge invariant perturbations 
can be written in terms of this Q we conclude that a^, are pure gauge 
terms which we can put equal to zero without loss of generality. 

Metric perturbations of Robertson-Walker space-times, which can be 
viewed as describing gravitational radiation, have also been studied by Bardeen 
|3] in an important paper. In this study the background space-time is taken 
to be a Robertson- Walker space-time with line-element 

ds^ = Q\T){-dT^ + '^gapdx'' dx^} . (7.1) 

Here the greek indices take values 1, 2, 3 and '^Qap is the metric tensor for a 
three-space of constant curvature. Comparing this to ()3.4|) we see that our 
background space-time also has this form if we take ^Qajs = {Po'^P, Po'^P, 1) 
and label the coordinates The method used in 

jS] involves separating the time dependent and spatial dependent parts of 
the perturbations. Now for us the important coordinates are u = {x — 
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T{t))/^/2, V = {x + T{t))/\/2 and there is no natural way to carry out this 
separation. Thus it is not possible to directly compare the results found 
here with those of [Sj. However there are some obvious similarities and 
differences between the results and we shall briefly comment on these now. 
One point of agreement is that gravitational radiation is described by tensor 
perturbations only. Specifically in our case gravitational waves are described 
by perturbations in the shear and anisotropic stress tensors. The perturbed 
space-time in |S| is given by 

ds^ = -n^dT"^ + gafsdx'' dx^ , (7.2) 

where 

9ap = n'i'gaf^ + 2H^^\t) qS(x^)] , (7.3) 

(2) 

and Qq,^ is a divergenceless trace-free tensor. This bears a strong resemblance 
to our perturbed space-time described by ()6.1|) where •jab given in ()6.2p is 
also divergenceless and trace- free. However it is clear from ()6.2|1 that, in 
effect, our small metric perturbations 'jab are expressible in the form of a 
2x2 matrix whereas {Q^^) is a 3 x 3 matrix. In addition ■jab satisfies the 
inhomogeneous wave equation ()6.12|) while Q^j satisfies the homogeneous 
wave equation [Hj 

g(2)"/3;7^^ + fc2g(2)a;3^Q^ (7.4) 

where is a constant. 
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A The Ricci Tensor Components 

In this section we give the Ricci tensor components (on the tetrad given by 
Eqs. ()5.8p ) for the metric defined by Eqs. ()5.8|) and ()5.9|) . In the calculation 
of the Ricci tensor components we use dVt/du = —dQ/dv and df/du = 
df /dv to simplify equations. Also for convenience we shall use subscripts 
y, z, u, V to indicate partial derivatives with respect to these variables and 
a prime to denote partial differentiation with respect to v. Following the 
Cartan method to find the Ricci tensor components we first find the non- 
zero Ricci rotation coefficients to be: 




(A.l) 
(A.2) 
(A.3) 
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^ f-'K(3z- n-' po - ^n-' r' {l-a)Ky 
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(A.4) 
(A.5) 
(A.6) 
(A.7) 
(A.8) 
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(A.9) 
(A.IO) 



We note that in this calculation we have discarded any terms which are 
second order or smaller in a, (3. Using these coefficients we obtain the Ricci 
tensor components: 



i?i2 = 2^-^^: {l3,-l3u)-2n-^ I3u.-2f-^ f n-\(3^ + (3,) , (A.18) 

+ 2fi-3fi'(a, - a„) - 2fi-2/"V'(a. + ««) - 2fi-'a„. 

+ Q-'f-'{2pt{p^'P)y, + Kzpl{pfP)y + Kypl{p^^P),} , (A.19) 

R22 = 2Q~"'Q'^ -2r' f" Q-^ + 2Q-^Q" -2Q-^f-^f'^ + Q-'f-^K 

- 2n-^n'{a, - au) + 2n-^f-'f'{a, + a„) + 2n^^auv 

+ Q'^f-^{p^{pQ^a)yy - PoiPo^a),, - KzplipQ^a), + Kypl{j)fa)y} 

+ ^-^f-\2pl{p-,^P)y, + Kzpl{p-,^P)y + Kypl{p~,^P),} . (A.20) 

B The Existence of Gauge Terms if , are 
non— zero 

In this Appendix we demonstrate that a^, appearing in ()5.6p are pure 
gauge terms. For clarity we shall consider only cases when / 7^ 0. When 
7^ 0, 7^ the line-element ()5.6p with 6, c given in ()5.7|) can be written 
in the form 

ds^ = 2n^dudv + pf fn'^{{l + a)dy + f3dz + Adu + PdvY 

+ p^'^fn'^{(3dy + {l-a)dz + Bdu + QdvY, (B.l) 




plr'^-^{{pYo^)yu + {pY ^)zu} , 
-plf-'n~^{{pfa),^-{pf^)y^} , 

4 ^'2 _ 2 ^]-3 _ 2 r V" fi"' , 

4^]-4^'2_2^]-3^]"-2rV"^^"' , 

41]-3Q"_2f]-4^]'2_2/-l/"^]-2 , 



(A.ll) 
(A.12) 
(A.13) 
(A.14) 
(A. 15) 
(A.16) 
(A.17) 



-plr'n-\{pYa),, - {pfP)y,} , 



where 



A 



al e" cosh (3 + a\e " sinh (3 , 
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B = a\e" sinh l3 + aj e"" cosh /? , 
P — ale"" cosh P + al e"" sinh l3 , 
Q = 4 sinh + g-" cosh P . (B.2) 

We find it convenient to work on the following tetrad: 

0^ = p-^fn{{l + a)dy + f3dz + Adu + Pdv}, 

= pQ^fn{(3dy + {l-a)dz + Bdu + Qdv}, 

6^ = ndu , 

= fldv . (B.3) 

As before our first step is to calculate the Ricci tensor components. In this 
case they are found to be: 

Ri3 = 2p^'f'n~\p^-A,)-p^'fn-'n'{Pu-A,) 

+ \po'm-\Pu - A,)u + pfr'^-\{pfa)yu + {pf(5).u} 

- \plr'^-^ (B^Kp^'y - ^BK^'y z - p^^B^^j 

- ]^plr'n-^ {pfA,,-pfKA,z + AKpf - ]^AK+'p^'y^^ , (B.4) 

R23 = 2p,'f'Q-\Q^-B,)-p,'fQ-^Q'{Q^-B,) 

+ ^p^'fn-\Q^ - B,)^ - Po'f-'^-'{{Po'a).u - {Po'/^)yu} 



Iplr'n-' {^AyKp^h - IaK%Sz - p.'Ay, 



- \plr'^-^{pfBy,-pfKByZ + BKpf-^-BK^p^'z^^ , (B.5) 

i?33 = 4n-'n''-2n-'n"-2Q-^f-'f"+pln-^{{p^^A)y^ + ip^^B),^} 

+ 2n-'f-'fpl{{p^'A)y + ip^'B),} , (B.6) 

i?44 = 4Q-^Q'^-2Q-^Q" -2Q-^f-'f" +p^^n-mp^^P)y, + {p^^Q),,} 

+ 2 fpliiPo'Q). + iPo' P)y} , (B.7) 

i?34 = 4Q-^Q" -2Q-*Q'^ -2Q-^ f-^ f" +plQ-^Q'{{p^^A)y + {po^B),} 

- pI n-' n'{{Po'P)y + {Po'Q)z} + \pI ^-\{pfA)y, + {pfB),,} 

+ ]^pI n~^{{pfP)y^ + + PI r' f ^-^{{pfA)y + {P,^B),} 
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+ plf-'f'^-H{Po'P)y + iPo'Q)z} , 



(B.8) 



Ri4 = plf-'n-'{{Po'a)y, + {p^'(3),,}-p^'fQ-'Q'{P^-A,) 

- 2p, 'f'Q-'{Pu - A) - ^Po 'f^-'iPu - 

- Iplr'n-' - P^Kp-' z + PKpf - \PK^P-,' y^] 

- \plr'^-^[QzKpfy-p^^Qy,-\QK^p^'yz^, (B.9) 

RiA = -plr^^-\{Po^a)zv - {Po^P)yv} - Po'm-'^'iQu - B,) 

- 2p^^f'n-'{Q^ - B,) - ^Po'fn-^iQu - B,), 

- \plr'n-^ [po^Qyy - QyKpf y + QKpf - ^QK^' z'^ 

- ^plr'Q-'I^PyKp^'z-p^'Py,-^PK%'yz^ , (B.IO) 

Ru = ^~^f~^{pt(Po'^(^)zz - pt(Po^(^)yy + KzpKpQ^a)^ - Kypl(p-^a)y} 

+ n-'r'{2ptipo'P)y. + Kzplip^'p)y + Kyplip^'p),} + n~'r'K 

+ ipln-^n' + plr'f'n-''){3p^'Ay - 2AKp^' y - 2BKpf z + p^'^B,} 

- ipln-'n' - plf-'f'n-'){3p-'Py - 2PKpf y - 2QKpf z + pfQ,} 

- ^Pq^Q-^K z{B^ + Qu) + (^Ay^ - ^A^Kpq^ y + Puy- ^PuKp^^ y^ 

- 2fi-2«„, - 2n-^f-^f'^ + 2Q-^n'2 - 2n-^r'f" + ^n-^'n" 

- 2n-''n'{a^ - a,) - 2n-^f-'f'{a, + a^) , (B.ll) 

R22 = ^~'^f~^{Po{Po^a)zz - Po(Po + KzpKpQ^a), - Kypl{p^'^a)y} 

+ n-'r'{2p'M'P)yz + Kyplipffi), + Kzpl{pf^)y} + ^-^f-^K 

+ {pln-^^ + plr'f'n'^){^pfB-, - 2BKpf z - 2AKpf y + p^Ay} 

- {pln-^^ - plr'f'^-'){^pfQz - '^QKpo' z - 2PKpf y + pfPy} 

- ^Po^fl-^Ky{A^ + Pu) + (^B,^ - ^B^Kp^^ z + Quz- ^QuKp^^ z^ 

+ 2Q-2a„„ - 2Q-V~ V + 2n-^n'^ - 2Q-V" V" + 

+ 2Q-^n'{au-a,) + 2n-T'f'{a, + au) , (B.12) 

Ri2 = n-^n' {A, + By) - n-^n' (p, + Qy) + 2n-^n' {p^ - - 2n-^Puv 



26 



- 2n-'f-'f'{(3^ + p,) + n-^f-'f'{A, + By) + n-^f-'f\Qy + p,) 

+ ^n-\A,, + By,) + ^n-\P,^ + Qyu) . (B.13) 

Here the subscripts y, z, u, v indicate partial differentiation with respect to 
these variables, differentiation with respect to v is denoted by a prime and 
K is the constant introduced in Eq. ()3.4j) . Using the above Ricci tensor 
components and Eqs. ()3.2p . ()3.3p . ()5.11|) it is easily checked that the first of 
the field equations ()2.11|1 is satisfied provided we choose A, B , P, Q to satisfy 
the Cauchy-Riemann equations 

(Po ' A), = {p-,^B)y , {pf A)y = -ip,' B), , (B.14) 

(Po 2 P). = iPo'Q)y , iPo' P)y = -(Po ' Q). ■ (B.15) 

Next with A,B,P,Q satisfying these equations we find from the remaining 
two equations in ()2.1H) that the conditions 

p^'AK = {Po'a)y^ + {Po'f3)zu, (B.16) 

p^^PK = (Po + (Po'/3).. , (B.17) 

p^'BK = -{Po'a),u + {Po'P)yu, (B.18) 

p^'QK = -{p~^a),, + ip^'P)y, , (B.19) 

Pu = A, , Qu = B, , (B.20) 

are sufficient to have Qa = and iTab = except for 

vTii = {^R-'Rpl + R-'f-'f'pl){3p^'Ay-2AKp,'y-2BKp^'z + p^'B,} 

- i^R-'Rpl - R-'f-'f'pl){^Po'Py - 2PKp,' y - 2QKpf z + pf Q.} 

- ^R-^ Kp^^ z{B, + g„) + R-^ (^Ay, - ^A, Kpq^ y + Puy - \Pu Kp^^ y^ 

+ R~'^f{Po{Po'^a)zz - Pl{Po'^a)yy + ^^^0(^0^")^ " ^l/Po(Po ^")?^} 

+ R-^f{2pl{pfP)y, + Kzpl{pfP)y + Kypl{pfP),} - 2R-' a„„ 

+ -^R-''R{a,-au)-2R-^r'f\a, + au), (B.21) 



7r,2 = {^R-^Rpl + R-^r'f'pl){^pfB,-2BKpfz-2AKpfy + p^^Ay} 
- {^R'^RPI - R'T'f'plH^Po'Qz - 2QKp^' z - 2PKpf y + p,^ Py} 
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- Kp^' z{A, + Pu) + R-^ (b,, - ^B, Kp^' z + Pu.- ig„ Kp^^ 

+ ^"V^{Po(Po ^")^^ - Po(Po + Kzpl{p^'^a), - Kypl{pQ^a)y} 
+ R~^f{2pl{pfP)y, + Kzpl{pfP)y + Kypl{pfP),} + 2R~^ 

+ + + + Qyu) + V2R-^R{p, - p^) - 2R-^ (3uv 

- 2R-^f-^f\(3u + ^) . (B.23) 

Following the procedure described in Section 5 we now construct vf = (tth — 
7r22 — 2i7ri2)/2. Using the conditions above to cancel terms we arrive at 

TT = ^R-^R{a,-au-i{(3,-(3u)}-2R~^r^f'{a, + au-i{(3. + (3u)] 

+ -^R~^R{Ay -B, + Q,-Py- i{A, + By-P,-Qy)}- 2R-^{auv - tPuv) 

+ R-^f-'f'{Ay -B,+Py-Q,- l{A, +By + P, + Qy)} 

~\~ 2 ^ {^y' Bzv ")~ Ryu Qzu ii,Azv -\- Byy -\- Pzu ~\~ Qyu)} • (B.24) 

We want to write vf in the form given in Eq. ()5.26|) for some analytic function 
Q. Before we try to do this we note that on account of the conditions Pu = A^, 
Qu = By we can write 

A = F^, P = Fy , (B.25) 
B = G^, Q = Gy, (B.26) 

for some functions F{y, z, u, v), G{y, z, u, v) which satisfy the Cauchy-Riemann 
equations 

ip^'F)y = -iPo'G)z , ip,'F)z = {p,^G)y . (B.27) 
Substituting these into Eq. ()B.24j) gives 

TT —-j=R R^F^y Guz ~l~ ^vz -^vy ^(^Ruz ~l~ ^(tj/ Rvz ^vy)} 

~^ R f f {Ruy G^MZ ~l~ -^vy ^(^Ruz ~l~ Cry^y -\- Fy^ ~\~ Gyy^^ 

+ R~'^{Fuyv — Guvz — i{Fuyz + Guvz)} + \/2R^^ R{ay — — i{(3y — 

- 2R-T'f'{ay + au- ^(/?, + - 2R-\auy - 2/?™) . (B.28) 
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In order to write vr in the required form we choose 



- /{"t- - + - . (B.29) 

Noting that D = V2dy we find that this G does satisfy Eq. ()5.26|) with 7f 
given by ()B.28|) and / 7^ provided a, P take the following form: 

a = ^{Fy - G,) + qiy, z,u) , = ^(^. + Gy) + r{y, z, u) . (B.30) 

Here q, r satisfy the Cauchy-Riemann equations 

(Po ^ Q)y = -(Po ^ r), , (po 2 g)^ = _ (^B.31) 

We remark that this is the first time we have made use of the fact that 
/ 7^ 0. If / =0 then a, (3 have a different form to that given in the last 
equation. Thus we emphasise that the analysis which follows does not apply 
if / =0. When a, /? are given by these equations it is straightforward to 
check using the various Cauchy-Riemann equations that Tiah is trace-free and 
the conditions (B.16)-(B.19) are identically satisfied. Substituting the above 
expressions for a, /3 into Eq. ()B.29j) yields 

Q = ^pf f{qu + iru) ■ (B.32) 

With k given by the first of Eqs. ()5.22|) it is trivial to show that Q satisfies the 
wave equation ()4.25|) . Now A, B, P, Q do not appear on the right-hand side 
of ()B.32j) and hence a^, do not contribute to Q. Thus since the perturbed 
shear and anisotropic stress can both be written in terms of Q we conclude 
that a^, are pure gauge terms. 
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